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Abstract. The relation between the globular cluster luminosity function 
(GCLF, diV/dlogL) and globular cluster mass function (GCMF, diV/dlogM) 
is considered. Due to low-mass star depletion, dissolving GCs have mass-to-light 
(M I L) ratios that are lower than expected from their metallicities. This has been 
shown to lead to an M/L ratio that increases with GC mass and luminosity. We 
model the GCLF and GCMF and show that the power law slopes inherently 
differ (1.0 versus 0.7, respectively) when accounting for the variability of M/L. 
The observed GCLF is found to be consistent with a Schechter-type initial cluster 
mass function and a mass-dependent mass-loss rate. 



1 Introduction 

Even though globular cluster systems (GCSs) are considered to have formed 
during mergers, the shape of the globular cluster luminosity function (GCLF, 
diV/d log L) differs fundamentally from that of young massive clusters (YMCs) 
in merging galaxies. The luminosity function of YMCs follows a power law with 
index —2 over the full mass range down to a few 100 Mq, while the GCLF 
peaks at ~ 10 5 L . This has been attributed to the ongoing tidal disruption 
of globular clusters ( GCs) and the re sulting destruction of low-mass, faint GCs 
( Fall &: Zhang 200ll ; IVesperini 200ll ). In studies considering this mechanism, 
the mass-loss rate of GCs is considered to be independent of their mass (corre- 
sponding to a disruption tim e oc M) and th e mass-to-light ratio (M/L) is 
assumed to be constant (e.g., iFall k, Zhang] l200ll ; I Jordan et al.ll2007l ). 

In order to trace back the merger history of galaxies using their globular 
cluster mass function (GCMF, d-ZV/dlogM), it is essential to obtain an accurate 
description for its evolution. Although the above studies reproduced the peaked 
shape of the GCMF, they did not account for two aspects of GC evolution: 



(1) The mass-loss rate of GCs does depend on cluster mass (iBaumgardtl bOQl; 
iBaumgardt & MakindT2003l : lLamers et alll2005l : lLarsenll2009l ) . correspond- 
ing to idis M* 7 with 7 ~ 0.7 (see Eq.Q]). This is due to the nonlinear scal- 
ing of the disruption time with the half-mass relaxatio n time (idis oc t^ 5 , 
Portegies Zwart et all Il998l ; iFukushige Heggiel |2000| ) . A lower value of 



7 means that the dissolution rate of low-mass clusters is slowed down with 
respect to massive clusters. 
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Figure 1. Left: The observed GCLF of Galactic GCs (|Harrislll996l ). Right: 
The inferred GCM F of Galactic GCs (histogram) using MjLy = 3 (as in 
iFall fc Zhangl 120011 ). Overplotted is our model MF with a mass-dependent 
mass-loss rate (solid line, see Section [2]) adopting a dissolution timescale to — 
1.6 Myr. The dash ed line shows the m odel for a cluster mass-independent 
mass-loss rate (as in lFall fc~Z hang 2 001[ ). Error bars are la Poissonian. 



(2) The M/L ratio o f GCs is not constant (iMandushev et al.ll l99ll;lBau mgardt & Makinol 
120031 : iKruiissenl I200SI : iKruiissen Mieskd 120091 ) due to the preferential 
ejection of faint, low-mass stars from dissolving GCs. Because low-mass 
GCs have on average lost a larger fraction of their initial masses, M/L 
increases with mass or luminosity. 

The effect of these aspects of GC evo lution on the relation between the GCLF 
and the GCMF has been considered bv lKruijssen Portegies Zwartl (|2009l ). The 
slope of the disrupt ion-dominated low-mass side of the GCMF is always equal 
to 7 (|Fall fc Zhandl200ll : lLamers et alJl200fil ). A mass-dependent mass-loss rate 
(7 = 0.7) therefore yields a different GCMF slope than a mass- independent 
mass-loss rate (7 = 1). The observed slope of the GCLF ~ 1, which is seem- 
ingly consistent with the latter if a constant M/L ratio is assumed (see Fig. [1]). 
However, the trend of increasing M/L ratio with luminosity implies that this is 
not necessar ily true because the slopes of t he G CLF and the GCMF fundamen- 
tally differ. iKruijssen &; Portegies Zwartl ( 20091 ) have shown that the observed 
GCLF is in fact consistent with a GCMF with a low-mass slope of ~ 0.7 when 

accounting for the trend of M/L r atio with luminosity. 

We revisit the calculations of IKruiissen &: Portegies Zwartl ( 20091 ) and in- 
clude a more detailed model for the evolution of the stellar mass function (SMF) 
within a star cluster. Previously, the low-mass star depletion was approximated 
by increasing the lower stellar mass limit of the SMF. In reality, the SMF evolves 
more gradually. This has been included in new models of the evolution of the 
SMF in dissolving star clusters (jKruiissen l2009h . which are used here. 



2 Star cluster models 



We u se the parameterised cluster model SPACE (jKruiissen & Lamersll2008l : lKruii sscn 
2009), which incorporates the effects of stellar evolution, stellar remnant pro- 



duction, dynamical dissolution and energy equipartition. The mass-loss rate due 
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Figure 2. Left: Modeled relation between M/Ly ratio and luminosity Ly 
(solid line). Diamonds mark the luminosities of observed Galactic GCs, and 
the dashed line denotes the constant M/Ly ratio that would be obtained if 
low-mass star depletion were neglected. Right: GCMF derived from the GCLF 
using a luminosity-dependent M/Ly ratio. The solid curve is the same as in 
Fig. [1] while the dotted curves represent models for (from bottom to top) 
log (^o/Myr) = log 1.6 -I- {—0.5, —0.25, 0.25}. Error bars are lcr Poissonian. 



to dissolution follows from i^is = t^M 1 as ( Lamer s et al. 2005): 



dis 



M 

^dis 



to 



(1) 



with to a constant that represents the rapidity of dissolution and depends on 
the environment. 

The evolution of the SMF is compute d by considering the ejection rate as 
a function of stellar mass ( Kruiis sen 2009!). The adopted method accounts for 
mass segregation and dissolution in a tidal field by using the timescale on which 
energy equipartition is reached for different stellar masses and by comparing the 
stellar velocities with the escape velocity. 



3 The inferred GCMF using realistic mass-to-light ratios 

The modeled relation between the M/Ly ratio and V-band luminosity Ly is 
shown in the left-hand panel of Fig. [2 Variations in di ssolution times cale and 



metallicity move this relation horizontally and vertically ( Kruijssenll2008l ) . When 
considering a real GC system, thi s induces scatter off the depicted relation. 
iKruiissen &; Portegies Zwartl ( 20091 ) considered this scatter in a detailed model 



of the Galactic GC system and concluded that a relation such as the one that 
is shown in Fig. [2] can still be used as a mean for the entire GC population. 

The GCMF that is derived from the observed Galactic GCLF using the 
relation between Ly and M/Ly is shown in the right-hand p anel of Fig . [H I t 
is compared to our modeled GCMF that has evolved from a ISchechter (| 19761 ) 



initial cluster mass function with index -2 and exponential cut-off mass M* = 
2.5 x 10 6 Mq. The modeled and observed GCMFs are fully consistent, implying 
that the mass- loss rate of GCs is indeed mass-dependent (t^ls ^ M 1 ', with the 
proportionality constant depending on the environment) and that there indeed 
exists a weak trend of M/L ratio with mass and luminosity. As such, the results 
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shown in Fig. [5] confirm the conclusions of Kruijssen &: Portegies Zwart ( 20091 ) 
with more detailed cluster models. 



4 Discussion 



Using a realistic model for the evolution of the SMF in GCs, we have shown that 
the slopes of the GCLF and the GCMF differ. As GCs dissolve, their masses 
decrease more rapidly than their luminosities due to the ejection of faint, low- 
mass stars. As a result, the slope of the GCLF is ~ 1 for luminosities below 
the peak luminosity, while the slope of the GCMF is ~ 0.7 for masses below 
the peak mass. This is consistent with a mass-loss rate that depends on cluster 

mass and on the environment 

The new cluster models ( Krui jssen 2009), in which we account for the chang- 



ing slope of the stellar mass func tion rather than shifting the lower st ellar mass 
limit, show that the results from iKruiissen fc Portegies Zwartl (120091 ) also hold 



when more detailed models are applied. This substantiates that care should be 
taken when comparing GCLFs and GCMFs. Because of the relatively low num- 
ber of GCs for which dynamical masses have been determined, we argue that 
the best method to compare them would be to model the GCLF rather than the 
GCMF, while accounting for the variability of the M/L ratio. This would al- 
low for a more accurate interpretation of the GCLF when studying correlations 
between its properties and the history of its galaxy. 

Acknowledgments. JMDK thanks East Tennessee State University for 
hosting an excellent conference and for financial support. The Leids Kerkhoven- 
Bosscha Fonds (LKBF) is acknowledged for supporting attendance to the con- 
ference. This research is supported by the Netherlands Advanced School for 
Astronomy (NOVA), the LKBF and the Netherlands Organisation for Scientic 
Research (NWO), grant numbers 021.001.038, 639.073.803, and 643.200.503. 



References 

Baumgardt, H. 2001, MNRAS, 325, 1323 
Baumgardt, H. & Makino, J. 2003, MNRAS, 340, 227 
Fall, S. M., & Zhang, Q. 2001, ApJ, 561, 751 
Fukushige, T., & Heggie, D. C. 2000, MNRAS, 318, 753 
Harris, W. E. 1996, AJ, 112, 1487 

Jordan, A., McLaughlin, D. E., Cote, P., Ferrarese, L., Peng, E. W., Mei, S., Villegas, 

D., Merritt, D., Tonry, J. L., & West, M. J. 2007, ApJS, 171, 101 
Kruijssen, J. M. D. 2008, A&A, 486, L21 

Kruijssen, J. M. D., & Lamers, H. J. G. L. M. 2008, A&A, 490, 151 
Kruijssen, J. M. D., & Mieske, S. 2009, A&A, 500, 785 
Kruijssen, J. M. D., & Portegies Zwart, S. F. 2009, ApJ, 698, L158 
Kruijssen, J. M. D. 2009, A&A, accepted, [arXiv: 0910 .4579] 

Lamers, H. J. G. L. M., Gieles, M., Bastian, N., Baumgardt, H., Kharchenko, N. V., & 

Portegies Zwart, S. 2005, A&A, 441, 117 
Larsen, S. S. 2009, A&A, 494, 539 

Mandushev, G., Staneva, A., & Spasova, N. 1991, A&A, 252, 94 

Portegies Zwart, S. F., Hut, P., Makino, J., & McMillan, S. L. W. 1998, A&A, 337, 363 
Schechter, P. 1976, ApJ, 203, 297 
Vesperini, E. 2001, MNRAS, 322, 247 



